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Abstract 

We show that a contribution to the total curvature perturbation may be due to the 
presence of flat directions in supersymmetric models. It is generated at the first os- 
cillation of the flat direction condensate when the latter relaxes to the minimum of 
its potential after the end of inflation. We also point out that, if the contribution to 
the total curvature perturbation from supersymmetric flat direction is the dominant 
one, then a significant level of non-Gaussianity in the cosmological perturbation is also 
naturally expected. 



One of the most successful predictions of the inflationary theory, the current paradigm for 
understanding the evolution of the early universe, is the redshifting of quantum fluctuations 
of the field driving inflation - the inflaton - beyond the Hubble radius, leading to an imprint 
on the background scalar (density) and tensor (gravitational waves) metric perturbations 
[1] that subsequently seeds structure formation. The latest confirmation of the inflationary 
paradigm has been recently provided by the five-year data from the Wilkinson Microwave 
Anisotropy Probe (WMAP) satellite [2] on the Cosmic Microwave Background (CMB) radi- 
ation anisotropy. 

Despite the simplicity of the inflationary concept, the mechanism by which cosmological 
curvature (adiabatic) perturbations are generated is not yet established. In the standard 
slow-roll inflationary scenario associated with a single inflaton field, density perturbations 
are due to quantum fluctuations of the inflaton itself. While this possibility is in agreement 
with present CMB data ^ l3j, it is not the only one. In the curvaton mechanism [1], the 
final curvature perturbation ( is produced from an initial isocurvature mode associated 
with quantum fluctuations of a light scalar (other than the inflaton), the curvaton, whose 
energy density is negligible during inflation and which decays much after the end of inflation 
(see also [5], where curvatons in supersymmetric theories are discussed). Further proposals 
invoke the inhomogeneity of the inflaton decay rate [6], inhomogeneous preheating [7] and 
the generation of the curvature perturbation at the end of inflation |8]. 

In this paper we point out that the generation of a flat spectrum for the curvature pertur- 
bation may be the general consequence of the presence of flat directions in supersymmetric 
theories. Let us briefly sketch how this can happen. In supersymmetric theories there exist 
many F- and D-term flat directions which are lifted because of the presence of the soft 
supersymmetry breaking terms in our vacuum, of possible non-renormalizable terms in the 
superpotential and of finite energy density terms in the potential proportional to the Hubble 
rate [9]. As a consequence, the field along the flat direction will acquire a large vacuum 
expectation value (VEV). After inflation, the condensate starts oscillating around the true 
minimum of the potential which resides at = 0. If the condensate passes close enough to 
the origin, the particles coupled to the condensate are efficiently created at the first passage. 
The produced particles become massive once the condensate continues its oscillation leaving 
the origin and may promptly decay into light relativistic stateil. This process allowing the 
generation of light states is called instant preheating [10]. Now, the key point is that the 
initial conditions for the flat direction when the oscillation starts may not be the same in sep- 



^Notice that this prompt decay of the x's effeciently removes them from the resonant band and no 
resonant preheating is expected from the oscillating flat direction. This is the reason why we consider 
particle production only at the first oscillation. 
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arate horizon volumes. This happens if the degree of freedom associated to the phase of the 
fiat direction is sufficiently light during inflation to be quantum mechanically excited. As a 
consequence, the condensate oscillates around the origin of its potential starting from slightly 
different values in different patches of the Universe. These different initial conditions give 
rise to fluctuations in the comoving number densities of the light relativistic states produced 
during the decay process after instant preheating and, ultimately, to CMB anisotropics. In 
this sense, supersymmetric flat direction provide a concrete and natural realization of the 
idea that the observed perturbations are associated to some underlying global symmetry 
which is slightly broken during inflation [7j. 

The generic potential for a supersymmetric flat direction cf) during inflation is given by 



V(^) = (ft^ - cHj) + (a^^|±^^" + h.c.) + |A|^ M^;!. (1) 

where c/, a/ and A are constants of fh and a^fn are the soft breaking mass terms of 

order the TeV scale, Hj is the Hubble rate during inflation, M is some large mass scale, pos- 
sibly of the order of the reduced Planck mass Mp ~ 2.4 x 10^^ GeV, and n is an integer larger 
than three. The if/-dependent terms are induced at the supergravity level by interactions 
between the flat direction and the inflaton fleld in the Kahler potential [TT] . 

For c/ > and Hi ^ m, the flat direction condensate acquires a VEV given by 



'I 



where /3 is a numerical constant which depends on a/, c/, and n. The phase 6i{'x,t) is very 
likely to undergo quantum fluctuations during inflation. Indeed, its mass is given by 

rrig = naj(3 cos{n9i + 9aj + 6\)Hj. (3) 

Here 9\ and 9aj are the phases of the coefficients A and a/. For the sake of simplicity we 
assume these phases to be zero in the following, our results also hold in the most natural case 
in which they do not vanish. For small enough values of the a/-parameter, me is smaller than 
the Hubble rate Hj during inflation. In fact, without any flne-tuning, the parameter a/ may 
be extremely small or even identically vanishing. In the extreme case in which the inflaton is 
a composite fleld, it will appear in the Kahler potential only through bilinear combinations 
and ai ~ Hj/Mp. In the case of D-term inflation [12] aj is exactly zero. The same occurs 
if we consider a flat direction which is lifted by a non-renormalizable superpotential term 
which contains a single fleld ip not in the flat direction and some number of flelds which 
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make up the flat direction [9], 

w = T^i^r'^. (4) 

For terms of tliis form, is non-zero along the flat direction, but W = Q along it. Examples 
of this type are represented by the direction ue which is lifted by = {X/M)uude, since 
= {X/M)uue is non-zero along the direction, and by the Que direction which is lifted 
by the n = 9 superpontial W = {X/M)QuQuQuHDee since F^^ = {X/M)QuQuQuee does 
not vanish [13]. If W = along the flat direction, no phase-dependent terms are induced. 
Alternatively, the superpotential may vanish along the flat direction because of a discrete 
i?-symmetry. In such a case, when W exactly vanishes, the potential during inflation has 
the form [9] 

= H}Mlf{\<t>\' /Ml) + H',Mlg{r/M;\ (5) 

and the typical initial value 0/ for the condensate may be 0{Mp), rather than Eq. ([2]), and 
a/ ~ {Hi/MpY'. All these considerations show that during inflation the phase 6{^,t) may 
be an effectively massless degree of freedom. 

In the post-inflationary era, the flat direction starts oscillating around = 0. At which 
frequency these oscillations take place depends crucially on the post-inflationary inflaton 
dynamics. If the inflaton is very weakly coupled, it will undergo a long period of oscillations 
around the minimum of its potential and eventually decay into radiation. If this happens 
when the Hubble rate is smaller than m, the flat direction will be anchored to the minimum 
of its potential till if ~ m/3 when it will start oscillating with a frequency of order of rh. 

On the other hand, the flat direction may oscillate around = with a much larger 
frequency. Indeed, the inflaton may release the energy stored in its potential very rapidly 
within a Hubble time. This is expected, for instance, if inflation ends through a rapid 
waterfall transition induced by a second field whose energy density dominates the energy 
density of the Universe in this phase [Tj. The potential of the flat direction will receive 
if-dependent corrections through the non-renormalizable couplings of the flat direction in 
the Kahler potential to the second field. In particular a mass squared term Cafif^|0p may 
be induced. If Caf is positive, the flat direction starts oscillating around the minimum of its 
potential with a frequency of the order of Hj if the waterfall transition is fascl. Similarly, a 
phase-dependent term [asiHjcj)^ / M^~^) may be generated through the coupling of the flat 



^ Notice that a large positive mass squared of the order of Hj may be induced even during the first 
stages of the radiation phase again by the non-renormahzable couphngs of the flat direction in the Kahler 
potential to the light relativistic fields </iiight [H]- Indeed, if there is a non-renormalizable coupling of the 
form |<^p|(/)iight|''/M^ and thermal effects generate the variance (^^ght) ~ then the flat direction will 
acquire a mass squared ^ T^/M"^ ~ Hj, where we have assumed that radiation is produced promptly after 
inflation and M is of the order of Mp. 



3 



direction to the second field driving the end of inflation. 

In fact, the flat direction may oscillate around = with a frequency much larger 
than Hj. This happens if the inflaton energy is released through a preheating stage |15j . 
Fluctuations of the scalar fields produced at the stage of preheating after inflation are so 
large that they can break supersymmetry much strongly than the inflation itself. These 
fluctuations may lead to symmetry restoration along flat directions of the effective potential 
so that the condensate moves at the very early stage of the evolution of the Universe, 

1 /2 

during the preheating era, with a frequency which can be as high as (m$Mp) ' , where m$ 
is the inflaton mass. Furthermore, a phase-dependent term can be induced by the large 
fluctuations of the scalar degrees of freedom generated at preheating, (m$0"/Mp [T5] . 

From now on, we will assume that the flat direction starts oscillating around the minimum 
of its potential right after the end of inflation with a frequency of the order of Hi (even though 
the reader should keep in mind that the frequency may be in fact larger), an initial amplitude 
0/ and a phase dependent term of order of {aoscH i(f)''^ / M^~^) . All the considerations made 
so far lead us to treat 0/, Oosc and n as basically free parameters. Indeed, is important to 
keep im mind that also the power n of the non-renormalizable terms lifting the flat direction 
might not be necessarily the same during inflation. As we described above, this happens if 
different non-renormalizable operators lift the flat directions during and after inflation. 

If the condensate passes sufficiently close to the origin, it can efficiently produce any 
state which is coupled to it. We generically call this state x (it niight be Higgses, squarks 
or sleptons) and suppose that it is coupled to the flat direction through the Lagrangian 
term h"^ |0|^ \xf- Its effective mass is therefore given by = + h"^ |0|^, where is 
the corresponding soft-breaking mass parameter. At the first passage through the origin, 
particle production takes place when adiabaticity is violated ^0\, \rh^\ /m'^ > 1. The fiat 
direction continues its classical oscillation, now giving a large mass to the produced quanta 
of order ~ /i^|0p. If they couple to some other light fermionic degrees of freedom, they 
can efficiently decay into these light states. 

The comoving number density of x particles produced during the first fiat direction 
oscillation 




where |0*| and |0=k| are the minimum distance and maximum speed of the trajectory with 
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Figure 1: The angular dependence of the minimum distance from the origin plotted in 
units of its initial value for Oosc = 0.003 and n = 4 (see the Appendix). The numerical 
evaluation (with /3 = 0.1 and A = 0.01) is represented by the continuous line, while the 
analytical estimate of Eq. ([71) is dashed. 



respect to the origin 

r (i±^) 

10*1 ^ 7raosc|0j| ^ „^ sin(ngj) (7) 

and 

10,1 ^ Hi 10,1(1 + (3^) Vl + 4(aosc/n) cos(n^/). (8) 

Eqs. ([7]) and ([8]) are obtained in a short Appendix and comparison of these approximations 
with the numerical evaluations are shown in Fig. [1] and Fig. [2l The non-adiabaticity 
condition |m^| /m^ 1 implies that |0|/(/i|0p) 1 or 




Let us now compute the curvature perturbation associated to the hght particles generated 
through the instant preheating stage. The presence of the phase-dependent term in the 
potential of the flat direction ([1]) during inflation violates the U{1) carried by the and gives 
n discrete minima for the phase of 0. The potential in the angular direction goes like cos{n6i) 
during inflation (the reader should remember that for simplicity we have set to zero the phases 
of the parameters aj and A). If the phase is a light degree of freedom during inflation, the 
field 6 does not sit at the minimum of its potential as it is quantum mechanically excited. 
It may acquire a random value, but constant over scales larger than the present horizon. 
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Figure 2: The angular dependence of the flat direction velocity at (f)^, (the minimum 
distance from the origin) plotted for Oosc = 0.003 and n = 4 in units of Hi\(f)i\ (see the 
Appendix) . The numerical evaluation (with (3 = 0.1 and A = 0.01) is represented by the 
continuous line, while the analytical estimate of Eq. ([8]) is dashed. 



Therefore, when the fiat direction starts oscillating, its initial condition varies from patch to 
patch. When light particles are generated through the instant preheating phenomenon, their 
abundance will not be uniformly distributed. On the contrary, isocurvature fluctuations of 
their number density are expected on superhorizon scales. These isocurvature fluctuations 
in the light fields are given by 

7r/i|0*p\ 6\<j)^\ 2nh\(l)^\^ 6\(f)^\ 



ght 



'flight 



6n. 



where, from Eqs. and (IHl) 

S\(f)* 



n 



cos{n6 



sm{n6i 
2a. 



sm{n6j 



1 + (4/ra)aosc cos{n6i) 



69. 



(10) 

(11) 
(12) 



These expressions lead to 

^^light 



sm{n9i 



a/1 + (4/n)aosc cosinOi) 1 4^/1 + {A/n)aosc cos{nei) 



a/1 + (4/n)aosc cos{n6i 



+ n cos{n9f) 



■ (13) 
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The fluctuation vanishes in the appropriate hmits, Oosc = 0: in these cases there is no depen- 
dence on the phase in the potential of the flat direction and therefore the condensate starts 
oscillating around the origin from the same initial condition throughout all the Universe. 
Supposing now that the inflaton field during inflation has generated a negligible amount of 
curvature perturbation during inflation, we can finally estimate the total curvature pertur- 
bation ( generated by the supersymmetric flat direction 

> Plight^ Plight Vlight 1 (^Plight 1 Plight ^^light /-..x 
Q = -■ Clight = -n^ : - ~Q = ~n ' y^^J 

Ptot Ptot Plight 'J Ptot -J Ptot ^Ught 

where 

Piight/Ptot ~ Px/Ptot ~ rrixnjptof (15) 

In Eq. the coefficient —1/3 assumes that the energy density of the Universe is dominated 
by the inflaton oscillations at the time of the generation of the curvature perturbation. If, 
on the contrary, relativistic degrees of freedom dominate at that instant of time, one should 
replace the factor —1/3 with —1/4. We stress that the curvature perturbation remains 
constant on superhorizon scales whatever the subsequent dynamics of the Universe is. Notice 
that may not exceed the energy stored in the flat direction, < p^. The fluctuations in 
9 take the form 

tt2 / U \ cos ne I -2t 

where Hj is the Hubble parameter at the time when the scale k exits the horizon at the 
value of the scale factor a^. The slow-roll parameter e = —H/H'^ accounts for the fact that 
during inflation the Hubble rate is slowly decreasing with time. The power spectrum of the 
phase therefore reads 

(U \ 2 / I. \ inajfS cos ne I -2e 
^) (^) • 

We can therefore estimate the maximum value of the curvature perturbation by using the 
maximum allowed value of the aosc-term and assuming ~ 

Cmax^ (^jpj cot{nei)6e, (18) 
correspoding to a maximum power spectrum of the curvature perturbation of the order of 
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where we have chosen M = Mp. Barring possible fine-tuning over the angle 6*/, we find for 
instance that for n = 4 and Hj ~ 10~'^Mp, the right amount of curvature perturbation V^^'^ ^ 
2.5 X 10^^ is obtained for A ~ 10^^, corresponding to Oosc ~ {XHi/Mpy^'^/{VTThsm{n6j) ~ 
10^^. Our findings indicate that a large (and possibly dominating) curvature perturbation 
may be obtained through the interplay of the infiationary and post-infiationary dynamics. If 
during inflation the field parametrizing the phase of the fiat direction has a mass smaller than 
the Hubble rate, then the particles generated through the phenomenon of instant preheating, 
when the fiat direction condensate oscillates around the minimum of its potential starting 
from slightly different initial conditions, will not be uniformly distributed on super-Hubble 
distances, thus leading to a nonvanishing curvature perturbation. This happens if particle 
production occurs. This requires that the flat direction passes sufficiently close to the origin, 
which imposes the bound ([9]) of the Oosc-term and, as a consequence, the upper bound (ITSl) 
on C- 

Let us close with some comments. The fact that we can choose non-renormalizable 
operators lifting the flat direction at the n = 4 level is relevant because there are, among 
all the flat directions lifted by n = 4 non-renormalizable operators, two of them, uude 
and QQQL, which carry no {B — L) number [9j (the ones which are lifted by n > 4 non- 
renormalizable operators all carry a non- vanishing {B — L)-number). This implies that the 
oscillations of the ude and QL directions may be associated to the right amount of curvature 
perturbation without giving rise to any (eventually large) baryon isocurvature perturbation 
associated to the baryon number generated by the flat directions [H]. Indeed, being uude 
and QQQL neutral under [B — L) any baryon asymmetry generated by their oscillations is 
promptly erased by the {B + L)-violating processes induced by sphalerons in the standard 
model. Notice that the operator QQQL may induce proton decay coming from the exchange 
of heavy Higgsinos unless we embed the MSSM within an SU(5) model (otherwise a flavour- 
independent Planck-suppressed QQQL with A ~ 10~^ will lead to rapid proton decay). In 
such a case, the scale M/A needs to be larger than about 10^'' GeV [19j . 

If (flSl) is the dominant component of the total curvature perturbation, the latter may 
have a sizeable non-Gaussian (NG) component which is easily found by expanding all the 
quantities obtained so far up to second-order in 59. The non-linear parameter /nl charac- 
terizing the level of NG [Hj becomes (3/nl/5) ~ — (Mp/0/)^(l/ cos^ 6i). Large values of NG 
are generically obtained unless the initial amplitude of the flat direction is close to Planckian 
values. Therefore, we conclude that, if the supersymmetric flat direction's dynamics produce 
a sizeable contribution to the total curvature perturbation, a large NG component in the 
CMB anisotropy is expected. 
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Finally, we point out that the curvature perturbation will be further diluited if, after 
its production, the energy density of the infiaton field dominates the energy density of the 
Universe for some period. We have to suppose therefore that the preheaing stage leads to 
full thermalization of the system. The possible problem related to the overproduction of 
gravitinos may be circumvented by a subseqyent late release of entropy [20] . 



Note Added 

When this paper was submitted to publication Ref. [18] was brought to our attention where 
the basic mechanism for the production of the curvature perturbation is identical to the one 
proposed in this paper. The decay mechanism of the fiat direction oscillations is however 
different and the possibility of generating a large NG was not discussed in Ref. jl8j . 
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Appendix 

In this short Appendix we will give an explicit derivation of the expressions of Eq. (JTj) and 
Eq. ([8]). Starting from the potential 



We rescale the fields in units of and time in units of the inverse frequency Hj^] the 
potential in units of becomes 

= + f + hA + (21) 



2' ' V n 

Since we work in the limit of small Oosc, we can write the real and imaginary part of the 
trajectories as 
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where 

~ cos(6'o)cost 



(0) ■ \ z ^^"^^ 

• ' ^ sm(6'o) cost, 



are the 0-th order (in aosc) solutions of the equations of motion (approximated for small (3) 
and t = (1 + P^)t. Solving the first order equations one finds 



- 2cos((n-l)^^o)/n(t)+0(aosc), 
-2sin((n-l)eo)/n(t) + 0(aosc 



.... ...w.. , ... . (24) 



where fn{t) is a complicated n-dependent function of time. As we shall see below, for the 
present purpose we are only interested in its value and the value of its derivative at t = 7r/2, 

,2 



\2) 2 r (1 + ^) ^ 



7r\ 1 , , 

25 

n 



2 r(l + f)^ ^"V2 

We can now analyse the approximated trajectories as they pass close to the origin, expanding 
the time-parameter around 7r/2 (the time when the 0-th order solution crosses the origin) as 
t = 71/2 + 6t. We then find that the minimum of the distance |0| from the origin lies at 

6i = -2ao3c/n (0 cos(neo) + OialJ. (26) 

Inserting this into the expressions of Eqs. (I22 i) -( l24l) . we can evaluate the minimum distance 
from the origin 



2aa 



/„,(!) I sin(n5o) + 0(aL), (27) 



which, when the proper dimensionful units are reintroduced, gives Eq.(j7]). Similarly one can 
evaluate the velocity taking the derivative of Eqs. (I22l) -( l24l) . and obtains Eq. (IHl). 
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